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We study WW and W Z production with luqq (I = e, /i) final states using data collected by the 
DO detector at the Fermilab Tevatron Collider corresponding to 4.3 fb~^ of integrated luminosity 
from pp collisions at y/s — 1.96 TeV. Assuming the ratio between the production cross sections 
a{WW) and a{WZ) as predicted by the standard model, we measure the total WV {V = W, Z) 
cross section to be a{WV) — 19.6 I3 q pb, and reject the background-only hypothesis at a level 
of 7.9 standard deviations. We also use &-jet discrimination to separate the WZ component from 
the dominant WW component. Simultaneously fitting WW and WZ contributions, we measure 
(t(WW) — 15.9 2 pb and a(WZ) — 3.3 j;* ^ pb, which is consistent with the standard model 
predictions. 

PACS numbers: 14.70. Fm, 14.70.Hp, 13.85.Ni, 13.85.Qk 



The study of the production of VV {V = W, Z) boson 
pairs provides an important test of the electroweak sector 
of the standard model (SM). In pp coUisions at -y/s = 1.96 
TeV, the next-to-leading order (NLO) SM cross sec- 
tions for these processes are a{WW) = 11.7 ± 0.8 pb, 
a{WZ) = 3.5±0.3 pb and a{ZZ) ^ 1.4±0.1 pb Q. Mea- 
suring a significant departure in cross section or devia- 



*with visitors from "Augustana College, Sioux Falls, SD, USA, 
''The University of Liverpool, Liverpool, UK, '^UPIITA-IPN, Mex- 
ico City, Mexico, '^SLAC, Menlo Park, CA, USA, '^University 
College London, London, UK, ■'^Centro de Investigacion en Com- 
putacion - IPN, Mexico City, Mexico, ^ECFM, Universidad Au- 
tonoma de Sinaloa, Culiacan, Mexico, and ''Universitat Bern, Bern, 
Switzerland. tDecoased. 



tions in the predicted kinematic distributions would indi- 
cate the presence of anomalous gauge boson couplings Q 
or new particles in extensions of the SM 0] . This anal- 
ysis also provides a proving ground for the advanced 
analysis techniques used in low mass Higgs searches 
The production of VV in pp collisions at the Fermilab 
Tevatron Collider has been observed in fully leptonic de- 
cay modes Q and more recently, in leptons-|-jets decay 
modes ^ , where the combined WW + WZ cross section 
was measured. 

In this Letter, we report observation of the production 
of a boson that decays leptonically in associated pro- 
duction with a second vector boson that decays hadron- 
ically {WV — ii^qq; £ = or 11^ , and v and q denote 
matter or anti-matter as appropriate) . The data used for 
this analysis correspond to 4.3 fb~^ of integrated lumi- 
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nosity collected between 2006 and 2009 by the DO de- 
tector at the Fermilab Tevatron Collider. The DO 
detector dijet mass resolution for W/Z decays of ~ 18% 
results in significant overlap ofW ^ qq and Z ^ qq dijet 
mass peaks. Therefore, we first consider WW and WZ 
simultaneously and measure the total WV cross section 
assuming the ratio of WW to WZ cross sections as pre- 
dicted by the SM. We then apply 6-jet identification to 
separate the WZ contribution, where the Z boson decays 
into bb pairs, from the dominant WW production. 

Candidate events in the electron channel are required 
to satisfy a single electron trigger or a trigger requiring 
electrons and jets, which results in a combined trigger ef- 
ficiency of (9813)% for the evqq event selection described 
below. A comprehensive suite of triggers in the muon 
channel, based on leptons, jets and their combination, 
achieves a trigger efficiency of (95 ±5)% for the fii^qq 
event selection. 

To select WV^ ivqq candidates, we require a single 
reconstructed electron (muon) with transverse momen- 
tum pt > 15 GeV (20 GeV) and pseudorapidity |ry| < 
1.1 (2.0) missing transverse energy > 20 GeV, and 
two or three jets reconstructed using a cone algorithm 
The jets must have pr > 20 GeV, I77I < 2.5, and at least 
two tracks within the jet cone Q originating from the pp 
interaction vertex. Lepton candidates must be spatially 
matched to a track that originates from the primary pp 
interaction vertex and they must be isolated from energy 
depositions in the calorimeter and other tracks in the 
central tracking detector. To reduce background from 
processes that do not contain W^ iv, we require that 
the transverse mass [13] is M^" (GeV) > 40 - 0.5.^t- In 
addition, we restrict M!^'^ < 200 GeV to suppress muon 
candidates with poorly measured momenta. 

Signal and most of the background processes are mod- 
eled with Monte Carlo (MC) simulation. The signal 
events are generated with pythia ^ using CTEQ6L1 
parton distribution functions (PDFs) [13] and include 
all SM decays. The fixed-order matrix element (FOME) 
generator ALPGEN [i3] with CTEQ6L1 PDF is used to 
generate M^-|-jets, Z+jets, and ti events. The FOME 
generator COMPHEP [iJ| is used to produce single top- 
quark MC samples with CTEQ6M PDF Both 
ALPGEN and COMPHEP are interfaced to pythia for par- 
ton showering and hadronization. The MC events un- 
dergo a GEANT-based [31 detector simulation and are 
reconstructed using the same algorithms as used for DO 
data. The effect of multiple pp interactions is included 
by overlaying data events from random beam crossings 
on simulated events. The next-to-NLO (NNLO) cross 
section is used to normalize the Z-|-jets (light and heavy- 
fiavor jets) (l6| . The approximate NNLO cross sec- 
tion [17| is used to normalize the tt samples, while the 
single top-quark MC samples are normalized to the ap- 
proximate next-to-NNLO cross section [l^. The nor- 
malization of the W^-|-jets MC sample (for all flavor con- 



tributions) is determined from data. Additional NLO 
heavy-flavor corrections are calculated with mcfm [1^ 
and applied to Z/VF-|-heavy-flavor jets MC samples. 

The multijet background in which a jet is misidenti- 
fied as a prompt lepton is determined from data. For 
the muon channel, the multijet background is modeled 
with data that fail the muon isolation requirements, but 
pass all other selections. For the electron channel, the 
multijet background is estimated using a data sample 
containing events that pass less restrictive electron qual- 
ity requirements. Both multijet samples are corrected for 
contributions from processes modeled by MC. The mul- 
tijet normalizations are determined from fits to the M^^ 
distributions and assigned uncertainties of 20%. 

To identify &-quark jets, in particular those originating 
from Z ^ bb decays, we use the DO neural network (NN) 
6-tagging algorithm [23|. The NN is trained to separate 
light-fiavor jets from heavy-flavor jets based on a combi- 
nation of variables sensitive to the presence of tracks and 
vertices displaced from the primary pp interaction vertex. 
The NN outputs for the two highest pr jets are then used 
as inputs to the flnal multivariate discriminant. We de- 
fine non-overlapping 0, 1, and 2-tag sub-channels based 
on whether neither, only one, or both of the two high- 
est Pt jets pass the least restrictive NN operating point, 
for which the &-jet identification efficiency and the light- 
fiavor jet misidentification rate are approximately 80% 
and 10%. Scale factors are applied to the MC events to 
account for any difference in efficiency or misidentifica- 
tion rate between data and simulation. 

The dominant background is VF-|-jets and therefore the 
modeling of this process in ALPGEN and the correspond- 
ing sources of uncertainties were studied in detail. Com- 
parison of ALPGEN with other generators {2l|] and with 
data shows discrepancies in jet 7y, dijet angular separation 
and the transverse momentum of the W boson candidate. 
Thus, data are used to correct these quantities in the 
ALPGEN W^+jets and Z-|-jets samples before ^-tagging is 
performed [22| . The possible bias in this procedure from 
the presence of the diboson signal in data is small, but is 
taken into account as a systematic uncertainty. 

As the diboson events are generated with a LO gener- 
ator, changes to the event kinematics and the acceptance 
due to a NLO and resummation effects are studied u sing 



events from the mc@nlo [23[ interfaced to HERWIG [24 1 
for parton showering and hadronization and using the 
CTEQ6M PDF set. Comparing kinematics at the gen- 
erator level after final state radiation, we parameterize a 
two-dimensional correction matrix in the pt of the dibo- 
son system and of the highest pr boson. After applying 
this correction to our pythia sample, we find good agree- 
ment with mc@nlo for all distributions studied. Half 
of the difference between the pythia and mcnlo pre- 
dictions is used as systematic uncertainty on the dibo- 
son production model, accounting for the possible effects 
of higher order corrections beyond NLO and of different 
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TABLE I: Number of events for signal and each background 
after the combined fit of WV using the RF output distribution 
(with total uncertainties determined from the fit) and the 
number of events observed in data. 





Electron channel 


Muon channel 


Diboson signal 


1725 ± 84 


1465 ± 67 


W/Z+light-RsLVOT jets 


37232 ± 1033 


33516 ± 709 


W/^+heavy-flavor jets 


5371 ± 608 


4854 ± 490 


tt and single top 


1746 ± 127 


1214 ± 86 


Multijet 


10630 ± 1007 


1982 ± 384 


Total predicted 


56704 ± 635 


43031 ± 531 


Data 


56698 


43044 
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FIG. 1: (color online) A comparison of the measured WV sig- 
nal (filled histogram) to background-subtracted data (points) 
in the RF output distribution (summed over electron and 
muon channels, and 0, 1, and 2-tag sub-channels), after the 
combined fit to data using the RF output distributions. Also 
shown is the ±1 standard deviation uncertainty on the back- 
ground prediction. The fit probability, P(x^)i is based on 
the residuals using data and MC statistical uncertainties. 



showering scenarios. 

The signal and the backgrounds are further separated 
using a multivariate classifier to combine information 
from several variables. This analysis uses a random forest 
(RF) classifier [Ulli, from which the output distribu- 
tion is used as a final variable to measure the produc- 
tion cross sections by performing a template fit. Fifteen 
well-modeled variables [27| that demonstrate a difference 
in probability density between signal and at least one 
of the backgrounds are used as inputs to the RF. The 
RF is trained using a fraction of each MC sample. The 
remainder of each MC sample, along with the multijet 
background samples, is then evaluated by the RF and 
used in the measurement. 

Depending on the source, we consider the effect of sys- 
tematic uncertainty on the normalization and/or on the 
shape of differential distributions for signal and back- 
grounds 27]. Systematic efi^ects on the differential distri- 
butions of the ALPGEN VF-l-jets and Z-t-jets MC events 
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FIG. 2: (color online) Results from the simultaneous fit of 
a{WW) and a{WZ) using the RF output distributions. The 
plot shows the best fit value with 68% and 95% confidence 
level (CL) regions and the NLO SM prediction. 



from changes of the renormalization and factorization 
scales and of the parameters used in the MLM parton-jet 
matching algorithm [28[ are also considered. Uncertain- 
ties onPDFs [29^, as well as uncertainties from object re- 
construction and identification, are evaluated for all MC 
samples. 

The total WV cross section is determined from a fit 
to the data of the signal and background RF output dis- 
tributions. This fit is performed simultaneously on the 
distributions in the electron and muon channels, and in 
the 0, 1, and 2-tag sub-channels, by minimizing a Poisson 
function with respect to variations in the systematic 
uncertainties [s^ . The magnitude of the systematic un- 
certainties is effectively constrained by the regions of the 
RF output distribution with low signal over background 
ratio. A Gaussian prior is used for each systematic uncer- 
tainty. The effects on separate samples or sub-channels 
due to the same uncertainty are assumed to be 100% cor- 
related. However, different uncertainties are assumed to 
be mutually independent. 

The fit simultaneously varies the signal and W^-l-jets 
contributions, thereby also determining the normaliza- 
tion factor for the W-|-jets MC sample. This obviates the 
need for using the predicted ALPGEN cross section, and 
provides a more rigorous approach that incorporates an 
unbiased uncertainty from Vl^-|-jets when extracting the 
signal cross section. The M^-|-jets normalization factor 
from the fit is consistent with the theoretical NNLO pre- 
diction |31i]. The yields for signal and each background 
are given in Table HI Though the total diboson yield 
includes a small contribution from ZZ — >■ ££qq events 
(1.5%), in which one of the charged leptons escapes de- 
tection, the cross sections presented here are corrected 
for this contribution assuming that the ratios between 
WW, WZ and ZZ cross sections are given by the SM. 
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FIG. 3: (color online) A comparison of the signal+background prediction to data in the RF output distribution (summed over 
electron and muon channels) for 0, 1, and 2-tag sub-channels after the combined fit to data using the RF output distribution 
(LP denotes light partons such as u, d, s or gluon, and HF denotes heavy- flavor such as cc or bb). The systematic uncertainty 
band is evaluated after the fit of the total WV cross section in the RF output distribution. 




FIG. 4: (color online) A comparison of the measured WW and WZ signals (filled histograms) to background-subtracted data 
(points) in the dijet mass distribution (summed over electron and muon channels) for 0, 1, and 2-tag sub-channels after the 
combined fit to data using the dijet mass distribution. Also shown is the ±1 standard deviation uncertainty on the background 
prediction. The fit probability, P(x^)i is based on the residuals using data and MC statistical uncertainties. 



The fit of the total WV cross section using the RF 
output distributions yields a{WV) = 19.6 tfxi P^: cor- 
responding to an observed (expected) significance of 7.9 
(5.9) standard deviations (s.d.). Figure [T] shows the 
background-subtracted RF output distribution summed 
over all sub-channels after the fit. As a cross check, we 
perform the measurement using the dijet mass distribu- 
tion in place of the full RF output distribution ^3] • This 
measurement yields a WV cross section of a{WV) = 
18.3 g pb, consistent with that obtained using the RF 
output distribution. 

The fit is then performed with the signal divided into 
the separate WW and WZ components, which are al- 
lowed to fioat independently. The result of this si- 
multaneous fit of a{WW) and a{WZ) using the RF 
output distributions is shown in Fig. [21 It yields 
a{WW) = 15.9 t\l (stat) t^l (syst) pb and a{WZ) = 
3.3 (stat) (syst) pb. The RF output distri- 

butions for the 0, 1, and 2-tag sub-channels from this 
fit are shown in Fig. [3l This measurement is also ver- 
ified fitting the dijet mass distribution, which yields 
cf{WW) = 13.3 l2.2 (stat) (syst) pb and (7{WZ) = 
5-4 i2.6 (stat) (syst) pb. Figure [4] shows plots for the 
background-subtracted dijet mass after the dijet mass fit. 

We also perform a fit in which we constrain the WW 



cross section to its SM prediction with a Gaussian prior 
equal to the theoretical uncertainty of 7% [l[ . The fit of 
the RF output distribution yields a WZ cross section of 
a{WZ) = 6.5 ± 0.9 (stat) ± 3.0 (syst) pb with observed 
(expected) significance of 2.2 (1.2) s.d., and the dijet 
mass fit yields (j{W Z) = 6.7 ± 1.0 (stat) ± 3.9 (syst) pb 
with observed (expected) significance of 1.7 (0.9) s.d. As 
expected, now that a{WW) is constrained to the SM pre- 
diction, the fit requires a higher rate for WZ in order to 
account for the excess of signal-like events. 

In summary, we have measured the cross section for 
total WV production to be aiWV) = 19.6 1^;^ pb 
{V — W or Z) with a significance of 7.9 s.d. above the 
background-only hypothesis. This result demonstrates 
the ability of the DO experiment to measure a dijet sig- 
nal in a background-dominated final state directly rele- 
vant to low mass Higgs searches. Furthermore, we have 
used 6-jet tagging to measure the contributions from 
WW and W Z and measured the cross sections for the 
separate processes to be a(WW) = 15.9 Ig ^ pb and 
(t{WZ) — 3.3 l pb. Although we cannot yet claim 
3 s.d. evidence of a WZ signal in the tiyjj final states, the 
extracted WV and WZ cross sections are in agreement 
with the SM prediction and their precise measurement 
represents an independent test to new physics which 
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could manifest itself differently in different final states. 
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APPENDIX 

INPUT VARIABLES TO THE RANDOM FOREST 
CLASSIFIER 

Here we define the fifteen variables used as inputs to 
the RF classifier. The observed distribution for each vari- 
able is shown in Figs. [SI [5] and [7] along with the predicted 
distribution after the fit of the total WV {V = W, Z) 
cross section using the RF output distribution. 

The RF input variables can be classified into three cat- 
egories: (i) 6-jet identification variables, (ii) kinematics 
of individual final state particles, and (iii) kinematics of 
multiple final state particles. Several variables are cal- 
culated using the four-momentum of the dijet system, 
which we define as the sum of the four-momenta of the 
two highest pt jets. We also reconstruct aW ^ candi- 
date W^'^, from the charged lepton and the I^t- The neu- 
trino from the VF— >■ iv decay is assigned the transverse 
momentum defined by I^t and a longitudinal momentum 
that is calculated assuming the mass of the £v system is 
80.4 GeV. Of the two possible solutions, we choose the 
real component that provides the smaller total invariant 
mass of all objects in the event. 

• b-jet Identification Variables: 

The NN 6-tagger has 12 operating points charac- 
terized by different purities. Each jet is assigned 
an integer &-tag value based on the highest purity 
operating point that it passes. 

1. Max &-tag Value: The greater 6-tag value of 
the two highest px jets. The neural network 
&-tagger has 12 operating points of increasing 
purity and the &-tag value corresponds to the 
highest operating point satisfied by the jet (or 
zero if the jet did not satisfy any of the oper- 
ating points). 

2. Min 6-tag Value: The lesser 6-tag value of the 
two highest px jets. 

These variables are shown in Fig. [SJ both in loga- 
rithmic and linear scales. 

• Kinematics of Individual Final State Parti- 
cles: 

1. pt{£)'- The Pt of the charged lepton. 

2. pT{}eti): The highest jet pr- 

3. pTQetj): The second highest jet pt- 

4. If^T'- The imbalance in transverse energy de- 
termined from the energy measured in each 
calorimeter cell and then corrected for recon- 
structed muons, jets, and electrons/photons. 

• Kinematics of Multiple Final State Particles: 



1. Mjj-. The invariant mass of the dijet system 
reconstructed from the two highest pT jets. 

2. M^" = ^2p^ (1 - cos(A0(^,^t))): The 
transverse W mass reconstructed from the 
charged lepton and the If^r- 

3. Ht — PT(jeti) +pT(jet2): The scalar sum of 
the two highest jet ptS. 

4. p-i(dijet,jeti)^ 

|pr(jeti) X prQeti -hjetj)]: The magni- 
tude of the leading jet transverse momentum 
perpendicular to the dijet system in the rest 
frame of the W^'^ candidate; 

5. p55'(dijet,jet2) = IpVOeta) x prOeti jetj)]: 
The magnitude of the second-leading jet 
transverse momentum perpendicular to the di- 
jet system in the laboratory frame. 

6. kf- = Afl(jet„jet,) p;"ffi^ : The angular 
separation between the two jets of highest pr, 
weighted by the ratio of the transverse mo- 
mentum of the second-leading jet and a scalar 
sum of the transverse momenta of the W^'^ 
constituents. 

7. cos(Z(dijet, jet^)): The cosine of the angle be- 
tween the momentum vectors of the dijet sys- 
tem and the highest pr jet in the laboratory 
frame. 

8. cos(Z(VF^'^,jeti))J'^': Cosine of the angle be- 
tween the momentum vectors of the leading 
jet and the W^^ candidate, evaluated in the 
rest frame of the dijet system. 

9. Centrality: The scalar sum of transverse mo- 
menta of the charged lepton and all jets in the 
event divided by the sum of their energies. 



SYSTEMATIC UNCERTAINTIES 

Table |IT] gives the size of the systematic uncertainties 
for Monte Carlo simulations and multijet estimates. We 
consider the effect of systematic uncertainties both on 
the normalization and on the shape of differential distri- 
butions for signal and backgrounds. Although Table 
lists uncertainties for the diboson and VF-|-jets simula- 
tion, these uncertainties are not used when measuring 
the diboson signal cross section, for which the diboson 
and M^-|-jets normalizations are free parameters. How- 
ever, the size of the uncertainty must be specified when 
estimating the significance and when we constrain the 
cross section for WW production to its SM prediction in 
the fit. 



9 



- D04.3fb ' 
r (a) 



10-* = 



10' = 



D0 4.3 fb ' 
(c) 



-I- Data 

■ wv 

I I W/Z+LP 
I I W/Z+HF 
^1 Top 

I I Multijets 
Uncert. 





8 10 12 

Max b-tag Value 

Data 

■ wv 

I I W/Z+LP 
I I W/Z+HF 
^1 Top 

I I Multijets 
t Uncert. 



10 12 
Max b-tag Value 



> 



10^^ 



D04.3fb' 
(b) 



-I- Data 

■ wv 

H W/Z+LP 

■ W/Z+HF 

■ Top 

I I Multijets 
m Uncert. 




c 1200- D04.3fb 

CD 




8 10 12 
Min b-tag Value 

Data 

■ wv 

n W/Z+LP 
H W/Z+HF 

■ 'Top 
I I Multijets 

^ Uncert. 



10 12 
Min b-tag Value 



FIG. 5: (color online) Distributions of the &-jet identification variables used as inputs to the RF classifier (first two of fifteen) for 
electron and muon channels combined, with logarithmic ((a) and (b)) and linear ((c) and (d)) scales. To better show the WW 
and WZ signals the lowest bin is cut off in the distributions with a linear scale. The signal and background predictions and the 
systematic uncertainty band are evaluated after the fit of the total WV cross section in the RF output distribution. Definitions 
for each variable are provided in the text (LP denotes light partons such as u, d, s or gluon, and HF denotes heavy-flavor such 
as cc or bb). 



Ax^ FOR WV MEASUREMENT 

The statistical significance of the diboson signal yield 
from the fit to the data is estimated via analysis of the 
Ax^ curve obtained by fitting the data to the sum of 
background and signal templates as a function of the sig- 
nal rate. The results of this analysis are given in Ta- 
ble [IlTl Figure [5] shows how the of the fit changes as a 
function of the signal cross section when using either the 
dijet mass or the RF output distribution to measure the 
total WV cross section. 
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FIG. 6: (color online) Distributions of the variables (next eight of fifteen) used as inputs to the RF classifier for electron and 
muon channels combined, and before &-tagging. The signal and background predictions and the systematic uncertainty band 
are evaluated after the fit of the total WV cross section in the RF output distribution. Definitions for each variable are provided 
in the text (LP denotes light partons such as u, d, s or gluon, and HF denotes heavy-flavor such as cc or 66). 
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FIG. 7: (color online) Distributions of the variables (remaining five of fifteen) used as inputs to the RF classifier for electron 
and muon channels combined, and before fe-tagging. The signal and background predictions and the systematic uncertainty 
band are evaluated after the fit of the total WV cross section in the RF output distribution. Definitions for each variable are 
provided in the text (LP denotes light partons such as u, d, s or gluon, and HF denotes heavy-flavor such as cc or 66). 
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TABLE II: The RMS amplitude (in percent) of each systematic uncertainty in the RF output distributions for the signal and 
background predictions. The RMS amplitude is defined as: \/Yl"=o Pi^l I X^ILo Vi'-, where pi is the predicted number of events 
in bin i, Ai is the percent change in bin i when the uncertainty is varied by 1 s.d., and n is the number of bins. In cases where 
the amplitude is different for different subchannels, the range of amplitudes is given. The rightmost column indicates whether 
the uncertainty only affects the normalization (N) or if it has also a differential dependence (D). 



Source of systematic 

uiiccrtcLintry 


Diboson signal 


VF+jets 


Z+jets 


Top Multijet 


Nature 


nilppi'TTin t'rio'O'PT' /Tl) pfirir'ipTir*^?' 


±5 


±5 


±5 


±5 


N 


IVEuon trigger /ID efficiency 


±5 


±5 


±5 


±5 


N 


Tpi" ir]pTi'l"ifipa'l"inTi 




±1-2 


±1-2 


± <l-2 




Jet energy sca-le 


±2-4 


±6-8 


±4-12 


±2-3 


D 


Tot onorcv rptjnliitinn 


±2-3 


±3-12 


±4-10 


±1-2 




Tpt vpvtpY pnn fi Tin ?i t i on 


±2-3 


±3-4 


±3-5 


±1-3 




Taggability correction 


± <1 


± <1 


± <1 


± <1 


D 


6-tae;e:ine: 


±1-5 


±1-4 


±1-5 


±8-10 


D 


Luminosity 


±6.1 


±6.1 


±6.1 


±6.1 


N 


Cross section 


±7 


±6.3 


±6.3 


±10 


N 


l/+heavy-flavor cross section 




±20 


±20 




N 


V+2 jets/V^+3 jets cross section 




±10 


±10 




N 


Multijet normalization 








±20 


N 


Multijet shape, electron channel 








± <1 


D 


Multijet shape, muon channel 








± <1 


D 


Diboson modeling 


±2-3 








D 


Parton distribution function 


±1 


±2 


±1-3 


±2-4 


D 


Unclustered Energy correction 


± <1 


± <1 


± <1 


± <1 


D 


ALPGEN jet r\ corrections 




± <1 


± <1 




D 


ALPGEN AR{jj) and pt{W) corrections 




± <1 


± <1 




D 


Re-weighting diboson bias 




± <1 


± <1 


± <1 


D 


Renormalization and factorization scales 




± <1 


± <1 




D 


Underlying event model 




± <1 


±1 




D 


ALPGEN parton-jet matching parameters 




± <1 


± <1 




D 



TABLE III: Results from fitting the total WV cross section and the uncertainties resulting from limited data statistics (stat), 
and possible systematic biases (syst). Also, the expected and observed Ax^ obtained by fitting the data with and without the 
specified signal process and the corresponding significance in number of standard deviations (s.d.) for a one-sided Gaussian 
integral. 





Measured a{WV) 


[Pb] 


Ax^ (significance) 

Expected 


Observed 


RF Output 


19.6 ± 1.4 (stat) +2.7 


(syst) 


35.8 (5.9 s.d.) 


63.5 (7.9 s.d.) 


Dijet Mass 


18.3 ± 1.5 (stat) +11 


(syst) 


22.1 (4.6 s.d.) 


33.0 (5.6 s.d.) 
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